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Abstract 

We study the global existence and space-time asymptotics of solutions 
for a class of nonlocal parabolic semilinear equations. Our models 
include the Nernst-Planck and the Debye-Hiickel drift-diffusion sys- 
tems as well as parabolic-elliptic systems of chemotaxis. In the case of 
a model of self-gravitating particles, we also give a result on the finite 
time blow up of solutions with localized and oscillating complex- valued 
initial data, using a method by S. Montgomery-Smith. 

1 Introduction 

In this paper we are concerned with semilinear parabolic systems of the form 

d t Uj = Auj + V • (XXfc=i c i-h,k u h {S7E d * u k )) , j = 1, . . . , m, (1) 
u(0)(x) = u (x). 

Here the unknown is the vector field u = (ui, . . . , u m ), defined on the whole 
space M. d (with m > 1 and d > 2), and Cj^,k £ -^°°(lR a! ), j,h,k = l,...,m, 
are given coefficients. Moreover, Ed denotes the fundamental solution of the 
Laplacian in M. d . 

Systems of the form ((TJ arise e.g. from plasma, semiconductors and elec- 
trolytes theories, biology (modelling of chemotaxis phenomena) and statisti- 
cal mechanics. The basic example for us is the model for gravitating particles: 
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in this case m — 1, d > 2, and the governing equations are usually written 

as 

d t u = Au + V- (uV<p), Aip = u. (2) 

Here u = u(x,t) is the density of the particles and (p is the self-consistent 
gravitational potential generated by u. Related systems appear also in the 
theory of chemotaxis, see e.g. jHj, |Zj, We do not require that u > in 
our study, which is, however, relevant in physical applications; we even admit 
complex- valued solutions. In this case the coefficients (cj^k) are constant and 
equal to 1. Another important example is provided by the Debye system, in 
which the first equation of (j2J is replaced with 

d t u = Au - V • (uVf). (3) 

A more general model, still belonging to the class (0), is the drift-diffusion 

system 

d t v = Av - (V • (vV4>)), 

d t w = Aw + (V- (wV(f))), (4) 
A(f> = v — w. 

In the theory by W. Nernst and M. Planck, v and w represent the density of 
positively and negatively charged particles, respectively. 

A lot is known about the existence and the nonexistence of real-valued 
solutions of these models, see e.g. 0, and the references therein. 

For instance, if d = 1, then the considered models have global in time so- 
lutions. If d > 2, the Debye system Q and more general (J1J) have global 
in time solutions and their asymptotics is described by suitable self-similar 
solutions, [Sj and [H]. These may be interpreted as a complete diffusion of 
charges to infinity due to repulsive interactions. 

On the other hand, models describing either chemotaxis or gravitational 
interaction in d > 2 dimensions feature concentration phenomena which may 
eventually lead to a collapse of solutions. These phenomena manifest by the 
formation of singularities of solutions like weak convergence either to Dirac 
point masses or to unbounded functions ~ \x\~ 2 - 

One purpose of this paper is to show that a different kind of finite time 
blow up can occur for solutions of © (and for a few other particular cases 
of Q). In particular, we will show that also nonpositive (in fact: complex- 
valued) and oscillating solutions can blow up. Our second purpose is to give 
a global existence result for "small" solutions of (0). Such result will provide 
us with some decay profiles in space-time of solutions. 

The global existence result for "well localized" solutions can be stated as 
follows (see section |21 for a more general, and more precise, statement). 
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Theorem 1.1 Let d > 3. There exists rj > such that if 

\u (x)\ < * (5) 
(1 + \x\) 2 

then there exists C > and a unique solution u of ^ such that, for all 
x ER d andt> 0, 

C C 
\u(x,t)\ < and \u(x,t)\ < j-^. (6) 

Moreover, if all the coefficients Cj^h,k{x) are constant in ~R d , then the smallness 
assumption (JHJ), can be replaced by the weaker, scale invariant, condition 

ess sup |x| 2 |mo(x)| < rj. (7) 



It would be possible to establish similar decay profiles in space-time for 
the solution, with a spatial decay rate larger than two. In this case the 
decay rate as t — > oo is also increased, up to one-half of the decay rate as 
x — > oo (or up to the rate d/2 if the space decay rate is larger than d). The 
exponent two is however the most interesting case, since it plays a special 
role in these models, for scaling reasons. For example, it corresponds to the 
expected decay rate of self-similar solutions, see e.g. It is also the decay 
of the well-known Chandrasekhar solution, u(x,t) = 2(d — 2)|x|~ 2 , which is 
a stationary solution for (J2J) for d > 3. 

We now state our result on the blow up of solutions. 

Theorem 1.2 There exists Uo G 5o(^ d ) (the space of functions belonging 
to the Schwartz class, with vanishing moments of all order), such that the 
corresponding solution u of (0) blows up in finite time: there exists t* > 
such that u(t*) B^(R d ) for all s el, 1 < p, q < oo. 

We will restate this theorem in a more precise way in section El Therein, 
we will also recall the definition of the Besov norm || • H^.^d-). Here we only 
observe that this theorem tells us, in particular, that ||-u(t)|| iP blows up for 
all 1 < p < oo. The proof of Theorem 11.21 consists in proving suitable lower 
bound estimates for the Fourier transform u. We shall derive such estimates 
using an idea of Montgomery-Smith [TU] . 

Our explosing solution u of Theorem 11.21 is in fact complex- valued since 
its Fourier transform enjoys some positivity and nonsymmetry properties. Of 
course, one can rewrite the scalar equation (0) for the real and imaginary 
part of u. This yields a blow up result for a real system of the form (fT|). 
which is formally close to (jlj). 

Notations. In chains of inequalities, all the constants will be denoted by 
C even if they vary from line to line. We will simply write J instead of J Rd . 
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2 Global existence for the general model 

The proof of Theorem 11.11 relies on size estimates of the kernel VE^. We 
have of course 

Our method applies also if we replace V-Ed with any kernel K such that K 
is a measurable function in M. d , and 

\K(x)\ < C\x\- d - 1+a , Ka<d. (8) 

In what follows, we will consider this more general situation. For the appli- 
cations that we have in mind, a = 2, and this explains the restriction d > 3 
in Theorem 11.11 The two-dimensional case is often a special case in these 
models, e.g. the Keller-Segel parabolic-elliptic model of chemotaxis. 

Since we would obtain the same bounds for all the components of u, in 
the remaing part of the section we can assume that u is scalar and Cj t h,k = c. 
Moreover, without loss of generality in the proof below, can assume that c 
is constant, essentially because multipling K with a L°° function does not 
affect ©. 

With this simplification, the models discussed above can be written in 
the following integral form 

u(t) = e tA u + [ G(t-s)*(u(K*u))(s)ds, (9) 
Jo 

where G behaves like a first order derivative of the Gaussian heat kernel, 
namely 

G{x,t) = t-( d+1)/2 V{x/Vi), with (10) 
e^G^-s^j =G(t-s). (11) 

To state our result in a precise way, we introduce a few useful spaces. 
For 9 > we define by the space of measurable functions / on M. d , such 
that (1 + | ■ \) 6 f G L°°(lR a! ). Let Eg the space of all measurable functions 
/ = f(x,t) in R d x R+, such that 

ess sup(l + \x\) e \f(x,t)\ < oo, (12) 



t>0 



ess sup(l +ty /2 \f(x,t)\ < oo, (13) 

t>0 



and 

feC{{0,oo),Lf). (14) 
The space £9 is equipped with its natural norm. 
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Theorem 2.1 Let 1 < a < d and K such that (JHJ) holds. Assume that 
u G . Then we can find r\ > such that if 

\\ u o\\l™<V, (15) 

then there exists a unique mild solution u of (JHJ), such that u G S a and 

u(t) — > Uq as t — > 0. Moreover, if we assume, in addition, that K is 
homogeneous of degree —d — 1 + a, then ()15|) can be replaced by the weaker, 
scale invariant condition: 

ess sup |x| a |Mo(x)| < 77. (16) 

x£R d 

The proof relies on the following simple lemma 

Lemma 2.1 Let K satisfy the assumption (JHJ)- 

1. IffeL™, thenK*feL™. 

2. If(f>e £ a , then K*<pe E x . 

Proof. Using the duality and the interpolation of Lorentz spaces, we get 
\\K*f\\»~ < ||A'|| £ -/(*fx-«).-||/|| x -/c-i),i < CII/II^LII/II^. 

Thus, 

\\K*f\\ L ~ <C||/|| L ». 
In particular, we may assume |x| > 1. Note that 

\K*f(x)\ <C [ \ X -y\- d -^\f(y)\dy = I 1 + I 2 + I 3 , 



where h = ...dy, I 2 = jM< M <f N • • • dy and I 3 = ...dy. 

One easily checks that these three integrals are bounded by C|x| _1 . The first 
part of the lemma follows. 

On the other hand, by the above inequality, 

wk * m\L~ < ^11011^^)11^)11^ < c i l + ^)" 1/2 ii^ii^- 

Combining this with the first part of the lemma, applied to 4>(t), yields the 
result. □ 

For u G £g, the nonlinear term u{K * u) belongs, by Lemma ETT1 to £q+\. 
Then it is natural to study the behavior of the linear operator 



L(w)(t)= I G(t - s) * w(s) ds (17) 

in such a space. 
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Lemma 2.2 Let 1 < a < d and w G S a +i- Then L(w) G £ a . 

Proof. We will use repeatedly the property 

\\G(t-s)\\ L ^=C(t-s)- 1 / 2 , 

which is a consequence of (jlOj) . A few estimates below bear some relations 
with those of Miyakawa |_9 i , yielding space-time decay results for the Navier- 
Stokes equations. We start observing that L{w) G L°°(R a! x R). Indeed, 

\\L{w){t)\\ L oo < [ ||G(*-s)|| L i|Ks)|| L <»<fc 
Jo 



< C\\w\\ £a+1 [ {t-sY^s^dsKCWwWe^. 
Jo 



Then we can assume in the following that \x\ > 1 and t > 1. 
We can write L{w) = Ii + I 2 , where, 

h= G(x - y)w(y,s) dyds 

Jo J\y\<\x\/2 

and 

ft 



I 2 = G(x - y)w(y,s) dyds. 

Jo J\y\>\x\/2 



Now, 



|/iOM)l < \A' d (t- s )- 1/2 (i + |y|)- a (i + s )- 1/2 ^rf s 

Jo J\y\<\x\/2 



< 



C\x\- d [ W a dy < C\x\ 



On the other hand, 

\h(x,t)\ < C\x\~ a I \\G(t - s)\\ L i s- 1 ? 2 ds < C\x\ 



Thus, \L(w)(x, t)\ < C(l + \x\) a \\w\\£ a+1 and, in particular, 

H-^( u ')llL-»((0,oo),L d /-.° c ) < C||w|| fQ+1 . 

To obtain a decay estimate as t — >• oo, we recall (fTT|) and write 
L{w)(t) = e tA/2 L(w)(t/2) + f G(t- s) * w ds = J x + J 2 . 

Jt/2 

By duality (we denote here by g t the Gaussian kernel), 

||Ji||loc < C\\g t /2\\ L d/(d- a ) > i\\L(w)(t/2)\\ L d/a,oo, 

< cr a/2 \\w\\ £ ... 

II II '-'CK + l 
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Moreover, 



\\J2(t)\\<Ct-^y*\\w\\ £a+1 f \\G{t-s)\\ L ,ds<Cr^\\w\\ £a+1 . 

Jt/2 

The decay estimates in space-time for L(w) then follow. The continuity with 
respect to t being straightforward, the proof of Lemma \2. 21 is finished. □ 

By Lemma 12. 2\ the bilinear operator 

B(u,v)= [ G(t-s)*(u{K*v)){s)ds (18) 
Jo 

is continuous from S a x S a to S a . Note that our last lemma also implies that 
\\u(t) — e* A Mo||L°° < C\ft, so that u(t) — > u a.e. and in the distributional 
sense. The existence (and the uniqueness) of the solution of (jSJ), under the 
assumption (|15jl now follows by a standard argument, i.e. the application of 
the contraction mapping theorem. 

In order to finish the proof of Theorem 12.11 it only remains to show that 
the smallness assumption (JT5J) can be relaxed, when the kernel K is a homo- 
geneous function. Consider the rescaling 

u x {x,t) = \ a u(\x,\ 2 t). (19) 

A direct computation shows that, if K is homogeneous of degree — d — 1 + a, 
and u is a solution of (jUJ), then u\ is a solution of (JOJ) as well. Now let r\ > 
be the constant obtained in the first part of Theorem 12.11 Assume that the 
datum uq is such that (fTB|) holds. Then we can choose a A > such that 

ess sup \ a {l + |x|) q |m (Ax)| < 77. 

We can apply the first part of Theorem 12. II to the initial datum A a w (A ■)■ ^ 
we denote by u the corresponding solution, we see that u^-i is the solution 
of starting from uq. Theorem 12. II is now established. □ 

Remark 2.1 With minor modifications of the decay exponents in the above 
proof, one sees that, for any finite T > the bilinear operator (j!8)l is bicon- 
tinuous in the space C([0, T], L^°), for all 9 > 0. The contraction mapping 
theorem guarantees that, if uq G 9 > (with arbitrary norm) and T > 
is small enough, then there exists a unique solution u G C((0, T], Lg°), such 
that u{t) — > uq in the weak sense; we will write u G C^QO, T], L^°) to express 
these properties. 
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3 Blow up for the model of gravitating par- 
ticles 



In this section we show that there exist solutions of (J2J), with initial data uq 
in the Schwartz class, and such that J x a ito(x) dx = for all a G N d , which 
blow up in finite time. Here we adopt a quite general definition of solution: 
we ask that the Fourier transform u(-,t), also denoted u t , satisfies for a.e. 
£ G R d and all t G [0, T], < T < oo, the integral equation 

u t (0 = e {s - tm2 M0 + T^p f e^ 2 it ■ (u s (0 * q|j5«.(0) da. W 

The definition of the Fourier transform for integrable functions that we adopt 
is = / u(x, t)e~~^' x dx. 

There are several ways to give a sense to the above integral and ensure 
the validity of (|2*UJl . An obvious way, is to consider the (local) solutions 
obtained in the setting of Remark 12. 1[ with 9 > d. But the above equality 
is true in more general settings. For example, it holds for the solutions 
u G C w ([0,T]: PM d ~ 2 ), (with < T < oo and d > 3) constructed in |U, 
where VM. a is the space of pseudomeasures 

VM a = {v G <S'(E d ) : v G L\ oc (R d ), \\v\\ VM * = esssup \£\ a \v(£)\ < oo}. 

As pointed out in [3], a distributional solution of the Cauchy problem for (J2J), 
does also satisfy (|2~Uj) . 

In this section we will consider initial data with nonnegative Fourier trans- 
form. Under this condition, one immediately checks that the iteration scheme 
yielding a solution in C w ([0, T], VM d ~ 2 ), or in C w ([0,T\,Lf), converges in 
the subset of functions u such that u(£,t) > 0, for all t G [0, T] and a.e. 
£ G M d . The crucial fact that will lead to the blow up is the following: 

Lemma 3.1 Let 

H 3 (u)(^t)= [ [ e^^^u^-V^s^dvds, j = l,...,d. 
Jo J \V\ 

Then, if0<u<v, suppu and suppv are contained in {£ G M. d : & > 0, i = 
1, . . . , d}, then < Hj(u) < Hj(v), and their supports are still contained in 
{£GM d :^>0,£=l,...,d}. 

This simple observation allows us to adapt to our situation the argument 
introduced by Montgomery-Smith for the "cheap" Navier-Stokes equations, 
see |H3|. 

Let us first recall the definition of the Besov norm II • II 00,00, a G R. 

^ 11 ni^j j 

We consider a function ip G 5(lR d ), such that ip > in R d , > 1 for 

I < |£| < 1, — for |£| < i or |£| > |. Then, for a distribution /, we 
define 

\\f\\^ = sup 2^ fe ||^(-2 fc )*/|| i0 o. (21) 

fcez 
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Remark 3.1 It is well-known that any Besov space Bp ,q (M. d ), as well as any 
Triebel-Lizorkin space F£ ,q (M. d ) (so in particular the L p -spaces, which are 
identified to -F°' 2 ), with s G R, 1 < p, q < oo, are embedded in Bt^ d ^ p '°° (M. d ) . 
It is then sufficient to show that u(t) blows up in the B^°° norm, for all a G R, 
to deduce that all Besov and Triebel-Lizorkin norms of u must blow up. To 
be more precise, L 1 is not a Triebel-Lizorkin space, but we will see that u t * 
becomes unbounded for a finite t*, hence ||w(£)||£i does blow up. 

A similar remark applies to pseudomeasure norms, since VA4 a is contin- 
uously embedded in B%~ (l ;x '■ m ''''' 
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Theorem 3.1 Let w G S(R d ), such that Wq is nonnegative and supported 
in the ball JBi/4(3ei/4), where e\ is the unit vector, and HwoIIl 1 = 1- Let 
A > 2 4//3 (27r) d and Uq = Awq (so in particular uq G iSo(R d )j. Assume also 
that u(-,t) is a tempered distribution such that for allt > 0, ut > and (120)1 
holds for a.e. £ G R d . Then, for all a G R, 

H-.OIUs 00 = °°' where t* = \og(2 1 ^). (22) 
Proof. Set to = 0, t k = log 2 (SjLi 2~ 2j j and w k = ■ We set also 

a*(f) = 2* k +^ k e-* H l t> _ tk (k G N), 
and claim that, for k — 0, 1, . . ., 

«t(0 > A 2fc « fe (t)^(0- (23) 
This is seen by induction. For k = the claim follows from Lemma f3. II 

ut(0 > Ae-W 2 w (0 > Ae-'woiO, t > 0. 
Now assume that (J23*Jl holds for k — 1. Set 

-E/c = {£ G R d : 2 fc ^ < 6 < |f | < 2 fc }, fc = 0, 1 . . . 

Note that w k = (2ir)^ d w k _i * w k _i. Thus, suppw^ C E k . 

But, for a.e. £ G E k , estimating from below by zero all the terms on the 
right hand side of fl2()L except for the first term obtained after computing 
the scalar product, we get 

Jt.JneE^ M 2 (^) d 
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In the second inequality we used our induction assumption. Now, for all 
t > t k , we have t-t k -i > 2" 2fe log2, so that 1 - e (*fc-i~*) 22fc > 1. This in turn 
implies f* i eC*-*) 22 * ds > 2~ 2k ~ 1 . Then, for all t > t k , we get 

u t (0 > A 2k 2 2k+7 ~ 2k e- 2kt , 

and ([SH)) follows. 

Moreover, ||wfc||n = (27r)~ d || 

Wk— 1 r l • 

Since 1 1 £y 1 1 x 1 — 1, by induction we 

get 

\\w k \\ L i = (2vr)- d ( 2fc - 1 ). 
Set t* = limfe^ootfc = log(2 1 / 3 ). We have ^(2~ k ■) > 1 in £ fc . Hence, 

^(2- fe 0^(0 > A 2k a k (t*)w k (0 > 0. 

Then 

)(a+d)k | 

•>afc I 



I^H^.oo > sup2 ia+ ^ fc |^(2 fc -)*^ 

fcgN 



= (27r)- d sup2 afc ||^(2- fc -)^IU 1 

fceN 

> su P A 2 ' : (27r)- d2 ' : 2 afc+2fc+7 - 2 ' : e- 2fer 

feGN 

= sup(fA2- 4 / 3 (27r)- d V 2 (" +2 ) fe+7 



OO. 



□ 



Remark 3.2 The same proof goes through for solutions of under the 
additional conditions Cj t h,k > for all j, h, k = 1, . . . , m, and, say, ci 7 h,k > 
c > 0. In this case, one can obtain a solution u = (ui, . . . , u m ), such that the 
first component blows up, starting from a datum uq such that «o > and 
with the first component satisfying the conditions of Theorem 13.11 

Remark 3.3 Analogous results can be obtained for space-periodic solutions 
of ((H), i.e. those defined on <i-dimensional torus. Instead of the Fourier 
transform, we will consider the Fourier coefficients u{£), (Gl 
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